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Since  1960,  tne  Defense  Research  Laboratory  has  been  working  on  problems 


00 


related  to  ref lec tion\and  scattering  of  sound  from  the  ocean  boundaries.  This 


work  has  been  sponsored  bv  the  United  States  Navy  Bureau  of  Ships,  now  the 

\CT>  T* 

>VBij  Naval  Ship  Systems  Command.  These  studies  have  included  both  analytical  and 
experimental  investigations,  with  the  experimental  work  being  done  in  a labo- 


cK 


ratory  tank  which  has  provisions  for  precise  control  of  the  parameters  of  the 


experiment.  Full-scale  experiments  at  sea  are  expensive  and  time  consuming, 
and  particular  aspects  of  reflection  theory  are  often  obscurred  by  sea  or 
bottom  conditions  that  cannot  be  controlled.  The  small-scale  experiments  in 
the  acoustic  tank  with  model  surfaces  provide  information  rapidly  and  accurately 
to  compare  with  mathematical  models.  These  mathematical  models  can  then  be  used 
to  predict  the  reflection  and  scattering  of  sound  by  the  ocean  boundaries  for 
full-scale  conditions. 


The  laboratory  tank,  asxillustrated  in  SLIDE  1 is  eight  feet  in  diameter  and 


seven  feet  deep.  A dual  transducer  positioning  system  with  an  accuracy  of 
± 0.003  in.  in  the  three  spatial  coordinates,  and  0.1  deg  in  azimuth  and  tilt 


Cl-  allows  the  geometry  of  the  scattering  situation  to  be  accurately  controlled.  The 


CJ>  model  being  studied  is  suspended  on  a tray  from  the  top  of  the  tank.  SLIDE  2 

illustrates  a simplified  block  diagram  of  the  equipment  used  in  the  experimental 
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that  standing  waves  are  not  set  up.  To  ensure  several  cycles  of  the  signal  in^py  gg  jgyy 


work.  Signal  measurements  in  the  tank  are  made  with  short  sinusoidal  pulses  so-^rpjjq  pjM 
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the  pulse,  the  carrier  frequency  of  the  pulse  is  usually  higher  than  signals 
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used  in  full-scale  experiments.  The  dimensions  of  the  model  are  scaled  with 
the  frequency  of  the  pulse  to  correspond  to  full-scale  conditions. 

Currently,  the  scattering  of  plane  waves  of  sound  from  a family  of 
rough  random  surfaces  is  being  studied  analytically  and  experimentally  at 
the  Defense  Research  Laboratory.  This  model  is  a good  representation  of  the 
sea  surface  and  many  areas  of  the  ocean  floor.  These  model  surfaces  are  the 
last  of  several  models  to  be  studied.  The  first  model  investigated  was  a 
plane  layer  of  fluid  sediment  overlying  a layer  of  sand.  This  model  is  repre- 
sentative of  some  abyssal  plain  areas.  The  next  model  consisted  of  a plane 
layer  of  sediment  overlying  a layer  of  limestone.  This  model  is  representa- 
tive of  other  abyssal  plain  areas.  The  analytical  methods  described  by 
Brekhovskikh  were  found  to  agree  well  with  the  experimental  results.  SLIDE  3 
illustrates  the  parameters  for  the  model  with  the  sediment  overlying  sand, 
and  SLIDE  4 indicates  the  theoretical  and  experimental  results  for  this  model. 
The  abscissa  is  the  grazing  angle  and  the  ordinate  represents  the  reflection 
coefficient.  The  agreement  of  the  theoretical  and  experimental  results  is 
obvious.  SLIDE  5 illustrates  the  parameters  for  the  model  of  sediment  overlying 
limestone,  and  SLIDE  6 indicates  the  theoretical  and  experimental  results  for 
this  model.  Again,  the  results  agree  quite  well.  A more  complete  description 
of  the  experiments  and  the  results  have  been  reported  by  Barnard,  Bardin,  and 
Hempkins^  in  the  Journal  of  the  Acoustical  Society  of  America. 


The  next  model  studied  was  a pressure  release  sinusoidal  surface,  which 
is  most  representative  of  the  sea  surface  but  also  corresponds  to  some  ocean 
bottom  areas  such  as  the  Blake  Plateau.  The  scattering  of  sound  by  a sinusoidal 
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surface  was  first  studied  by  Rayleigh.  He  obtained  amplitude  coefficients  for 
normal  incidence  in  terms  of  an  infinite  set  of  simultaneous  equations. 

Rayleigh  was  able  to  solve  this  set  only  with  the  restriction  that  the  incident 
wavelength  be  large  compared  to  the  sinusoidal  surface  amplitude.  In  the  past 
decade,  several  papers  have  appeared  in  the  literature  concerning  the  scat- 
tering of  sound  by  sinusoidal  surfaces.  Most  of  the  theories  advanced  thus 
far  contain  assumptions  that  in  general  provide  inadequate  descriptions  of  the 
reflection  of  sound  from  a simple  sinusoidal  pressure  release  surface  unless 

severe  restrictions  are  placed  on  the  geometrical  parameters  of  the  sinusoid 

2 

and  the  frequency  involved.  In  a recent  paper,  Uretsky  described  an  exact 

method  of  solution.  He  made  no  simplifying  assumptions  concerning  the  frequency 

of  the  incident  radiation  or  the  geometrical  parameters  of  the  sinusoidal 

surface.  Techniques  for  numerical  computation  were  not  stressed  by  Uretsky 

3 

but  have  since  been  reported  by  Barnard  et  al.  The  sinusoidal  surface  studied 
experimentally  at  the  Defense  Research  Laboratory  was  made  by  cutting  a sinusoidal 
variation  in  a piece  of  styrofoam  about  120  cm. on  each  side  and  about  7 cm.  thick. 
This  cutting  was  done  with  a special  tool  on  a milling  machine  and  produced  a 
sinusoidal  surface  with  a peak  to  peak  amplitude  of  3.0  cm.  and  a surface  wave- 
length of  4.5  cm.  SLIDE  7 is  experimental  data  illustrating  the  variation  of 
the  specular  reflection  coefficient  with  grazing  angle  when  the  carrier  frequency 
of  the  pulse  is  100  kHz.  At  100  kHz,  the  ratio  of  the  incident  sound  wavelength 
to  the  peak  to  peak  surface  amplitude  is  about  0.5, and  the  ratio  of  the  incident 
sound  wavelength  to  the  wavelength  of  the  surface  is  about  0.34.  For  the  data 
shown  in  SLIDE  8,  the  frequency  of  the  incident  pulse  was  400  kHz  where  the  ratio 
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of  the  wavelength  of  the  incident  sound  to  the  surface  amplitude  is  0.12,  and 
the  ratio  of  the  incident  sound  wavelength  to  the  wavelength  of  the  surface 
is  0.09.  SLIDE  9 illustrates  the  comparison  of  theoretical  and  experimental 
results  for  the  situation  where  the  incident  grazing  angle  is  fixed  at  45  deg 
and  the  reflected  grazing  angle  is  varied  from  about  30  deg  to  120  deg  with 
the  carrier  frequency  of  the  pulse  being  100  kHz.  It  is  evident  that  the 
scattering  from  this  surface  is  similar  to  the  scattering  one  observes  from 
an  optical  grating.  Corrections  have  been  made  for  the  directionality  of  the 
hydrophone  but  not  for  the  projector  in  the  theoretical  curves.  Corrections 
for  the  directionality  of  the  projector  should  improve  the  agreement  between 
the  two  curves.  SLIDE  10  illustrates  the  specular  reflection  coefficient  for 
a grazing  angle  of  45  deg  as  the  carrier  frequency  of  the  pulse  is  varied 
from  70  kHz  to  260  kHz. 

The  generally  satisfying  agreement  between  the  theoretical  predictions 
and  experimental  results  in  the  tank  work  and  other  comparisons  between  the 
results  of  the  model  tank  work  and  full-scale  results  at  sea  demonstrate  the 
usefulness  of  the  tank  work.  This  agreement  between  theory  and  experiment  for 
models  having  easily  defined  parameters  was  the  basis  for  continuing  the  model 
scattering  work  for  models  whose  parameters  were  not  so  easily  defined.  The 
remainder  of  this  paper  is  a description  of  such  models  and  the  theoretical 
and  experimental  results  obtained. 

The  sea  bottom  and  the  surface  do  not  in  general  have  the  regularity  of 
the  sinusoidal  surface  so  more  irregular  surfaces  were  considered.  The  model 
chosen  for  further  study  was  a rough  random  surface. 


A survey  of  contour  maps  of  the  ocean  bottom  indicated  none  were  sufficiently 
detailed  for  a model  so  aeromagnetic  maps  of  the  Northwest  Canadian  Shield  were 
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interpreted  as  topological  maps  and  yielded  surfaces  suitable  for  modeling. 

Four  of  the  maps,  each  40.6  cm  x 40.6  cm,  scaled  1.58  cm  to  the  kilometer, 


in  the  form  of  a square  were  used  as  the  pattern  for  the  model.  The  aeromagnetic 

.4 


maps  were  the  subject  of  a previous  study  by  Horton  et  al  which  determined 
their  statistical  properties.  This  earlier  work  concerned  applying  the  sta- 
tistical considerations  to  the  study  of  geological  studies. 

The  model  is  85  cm  in  length  along  an  edge  representing  a distance  of  51.5  km 
on  a map.  The  height  dimension  was  converted  from  the  aeromagnetic  dimension, 
the  gamma,  using  the  scale  10  y = 0.0795  cm.  The  contours  on  the  map  were 
located  on  the  model  with  an  accuracy  of  0.19  mm  in  each  of  the  three  coor- 
dinate variables.  This  is  an  accuracy  of  approximately  one  hundredth  of  a 
wavelength  at  100  kHz.  The  total  relief  of  the  model  from  the  highest  to  the 
lowest  points  is  slightly  less  than  four  wavelengths.  SLIDE  11  is  a photograph 
of  the  model  surface  used  in  the  studies.  In  order  to  avoid  difficulties  from 
the  edges  of  the  surfaces,  the  insonified  region  was  restricted  to  the  central 
fourth  of  the  surface  which  is  outlined.  The  elliptical  curve  shows  the  location 
of  the  half-power  level  of  the  incident  beam  when  the  incident  grazing  angle  is 
45  deg.  The  small  black  square  is  one  wavelength  on  each  side  at  100  kHz  so 
that  one  can  judge  the  size  of  the  hills  and  valleys. 

The  specific  formulas  for  the  scattering  coefficients  which  will  bt  used 
in  this  paper  to  compare  with  experimental  data  are  based  on  the  assumption  that 
the  statistics  of  the  rough  surface  are  isotropic.  In  order  to  test  this  assump- 
tion, the  autocovariance  function  of  the  surface  was  computed  and  the  values  are 
contoured  on  SLIDE  12.  In  order  to  make  this  computation,  a square  grid  of  sample 
points  % in.  apart  was  drawn  on  the  central  one-fourth  of  the  four  maps.  This 
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gave  1089  values  of  the  elevation.  The  autocovariance  function  was  computed  at 
intervals  of  \ in.  in  each  of  the  two  coordinates.  The  formulas  used  in  this 
calculation  have  been  given  elsewhere. 

The  contours  in  SLIDE  12  are  essentially  circular  for  small  distances 
while  they  become  elongated  ovals  for  larger  separations.  Fortunately,  the 
scattering  at  high  frequencies  is  determined  by  the  curvature  of  the  hills  and 
valleys  so  it  is  only  the  shape  of  the  contours  near  the  origin  that  is  important. 
Consequently,  one  may  assume  without  hesitation  that  the  surface  is  isotropic. 

Theoretical  Formulas  for  the  Model 

The  literature  on  the  reflection  of  waves  from  rough  surfaces  is  quite  exten- 
sive so  it  would  be  impossible  to  provide  a useful  guide  within  the  limits  of 
this  paper.  Fortunately,  a recent  book  by  Beckmann  and  Spizzichino^  provides 
an  excellent  survey  of  the  literature.  Of  the  many  theoretical  solutions 
available  in  the  literature,  it  was  decided  that  the  paper  by  Eckart^  offered 
the  greatest  promise.  His  solution  is  easily  evaluated,  and  it  does  not  require 
elaborate  mathematical  analysis  such  as,  for  example,  the  solution  of  a set  of 
infinite  equations.  Eckart's  assumptions  related  to  the  statistics  of  the  surface 
are  that  the  surface  heights  have  a guassian  distribution,  which  is  valid  for  this 
model  and  that  the  autocovariance  function  of  surface  heights  is  of  gaussian  form. 
The  autocovariance  functions  in  the  x and  y directions  for  this  random  surface 
model  are  illustrated  in  SLIDE  13.  As  is  evident,  the  functions  are  more  nearly 
exponential.  A straight  line  approximation  is  illustrated  through  the  first  four 
points  in  the  x direction. 

Although  the  statistical  properties  of  the  present  model  surface  do  not 
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satisfy  exactly  the  assumptions  on  which  Eckart's  theory  is  based,  the  agreement 
is  close  enough  to  justify  a comparison  between  experiment  and  theory.  Actually 
it  is  not  likely  that  the  departure  of  the  surface  statistics  from  gaussian  sta- 
tistics will  affect  the  theoretical  calculations  seriously  until  one  attempts  to 
calculate  standard  deviation  fora,  the  scattering  coefficient.  Another  reason 
for  the  appeal  of  his  solution  to  the  authors  is  that  it  is  based  on  Helmholtz's 

formula  which  has  provided  very  practical  formulas  for  other  applications  at 
7 8 

this  laboratory  * . Consequently,  it  was  decided  to  use  and,  if  necessary,  modify 

Eckart's  solution  unless  disagreement  with  the  experiment  required  that  it  be 

abandoned.  Fortunately,  the  agreement  is  sufficiently  good  that  there  is  no 

9 

incentive  to  try  other  theories.  A paper  by  Horton  and  Muir  contains  several 
explicit  formulas  derived  from  Eckart's  work,  and  a companion  paper^  provides 
more  detail  of  the  work  described  in  this  paper. 

The  Horton  and  Muir  paper  concerns  the  application  of  Eckart's  scattering 
theory  to  surfaces  with  various  statistical  properties.  Scattering  coefficients 
are  developed  for  both  the  low  and  the  high  frequency  cases  for  isotropically 
rough  surfaces  with  exponential,  gaussian,  and  sinusoidal  autocovariance  functions. 
The  specific  assumptions  of  interest  for  the  model  under  consideration  are  that 
surface  heights  have  a gaussian  distribution,  the  autocovariance  function  of 
surface  heights  is  exponential,  and  the  wavelength  of  the  incident  sound  is 
small  compared  to  the  maximum  relief  of  the  surface.  The  scattering  formula 
obtained  by  Horton  and  Muir  for  these  assumptions  with  a linear  approximation 
for  the  autocovariance  function  is  shown  as  Eq.  1 in  SLIDE  14, 
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where  the  parameters  of  the  equation  are  also  shown. 


ctq  = l/8n(k2hLc2)2 


k = radiation  wavenumber 


h = mean  square  surface  relief 


2 2 2 
w = (kL)  (a  +b  ) 


a,b,c  = sums  of  the  direction  cosines  of  the  source 
and  receiver 


L = autocovariance  decay  function  distance 
In  order  to  examine  more  easily  the  predictions  of  this  solution,  it  is 
desirable  to  introduce  the  actual  variables  used  in  the  experimental  studies. 


1 


These  variables  are  illustrated  in  SLIDE  15.  The  source  P was  always  held  at  a 
fixed  position  in  the  (x,  ) plane  while  the  receiver  R was  moved  by  varying  one 
of  the  two  angles,  0^  and  4^.  In  the  majority  of  measurements,  9^  was  held  con- 
stant equal  to  zero.  In  terms  of  these  variables,  some  of  the  direction  cosines 
are  shown  on  this  slide.  When  the  formula  for  a is  expanded  in  these  variables 
with  the  further  specification  that  0^=0,  the  expression  for  the  scattering 
coefficient  becomes  Eq.  2 on  SLIDE  16  where  G is  a dimensionless  constant  defined 
by  Eq.  3.  The  parameter  G characterizes  the  relief  and  the  coherence  of  the 
surface  in  terms  of  the  wavelength. 

The  subscript  T has  been  applied  to  a to  distinguish  this  theoretical  value 
from  the  experimental  value  a . It  will  be  shown  on  a subsequent  slide  that  this 

Ct 

formula  does  not  agree  very  well  with  the  experimental  observations.  In  particular, 


does  not  approach  zero  as  4^  — 0 or  tt,  and  the  maximum  in  a^,  considered  as  a 
function  of  0 does  not  occur  at  = 4 as  the  experimental  data  suggests. 
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This  serious  disagreement  between  theory  and  experiment  led  to  a renewed 
examination  of  the  assumptions  underlying  Eckart's  theory.  Mr.  Mitchell,  who 
performed  much  of  this  work  for  his  Master's  degree  thesis  noticed  that  sig- 
nificantly better  agreement  between  theory  and  experiment  would  be  obtained  if 
one  changed  the  boundary  condition  on  ap^/dn  in  the  Helmholtz  integral.  Eckart 
suggested  that  one  could  use  Eq.  4 on  SLIDE  17  on  the  pressure  release  surface, 
but  he  pointed  out  that  "in  the  deep  shadows  it  would  be  more  reasonable  to 
suppose"  the  condition  of  Eq.  5.  This  latter  condition  implies  that  the  pres- 
sure field  vanishes  at  the  surface  in  the  shadowed  regions.  On  the  surface  of 
strong  relief,  like  the  present  model,  both  conditions  exist  on  different 
portions  of  the  surface  so  that  one  is  in  closer  agreement  with  the  facts  if 
one  assumes  the  condition  of  Eq.  6 over  the  entire  surface. 

When  the  analysis  is  modified  by  this  boundary  condition  on  the  scattered 
wave,  P^,  the  only  change  in  the  expression  obtained  by  Horton  and  Muir  for  the 
general  integral  equation  for  the  scattered  pressure  field,  shown  as  Eq.  7 in 
SLIDE  17, 


._  ...  ikr10,  |’f_  -ik(ax+by+cC)  , 

4np^  = ikc(e  /r^  P e dx  dy 

is  that  the  factor  C should  be  changed  to  C where  C is  the  direction  cosine  of 

A A 

the  source  and  C is  the  sum  of  the  direction  cosines  of  the  source  and  the  receiver. 
When  this  change  is  followed  through  to  the  equation  for  the  scattering  coefficient, 
it  becomes  Eq.  8 in  SLIDE  17, 
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In  terms  of  the  variables  shown  in  the  previous  slide,  the  expression  for  the 
modified  scattering  coefficient  becomes  Eq.  9.  The  subscript  M is  added  to  a to 
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indicate  that  this  formula  is  obtained  from  a modification  of  Eckart's  theory. 

The  general  formula  without  the  restriction  0^  = 0 is  given  in  Eq.  10  where 
a,  b,  c,  and  C were  defined  on  a previous  slide.  The  organization  of  this 

A 

equation  has  been  changed  slightly  in  order  to  separate  the  angular  parameters 
from  the  parameter  G that  characterizes  the  linear  dimensions. 

Experimental  Procedures 

Only  a brief  account  of  the  experimental  procedures  will  be  given  here. 

If  further  information  is  desired  on  experimental  equipment  or  procedures, 
then  the  authors  would  be  pleased  to  provide  this  after  this  session.  The 
experiments  were  performed  with  the  model  surface  submerged  in  a water  filled, 
redwood  tank  24  ft  in  diameter  and  20  ft  high.  Since  it  was  easier  to  move  the 
hydrophone  in  a horizontal  plane,  the  model  surface  was  fixed  in  a vertical 
plane  on  a weighted  mounting  frame.  SLIDE  18  shows  the  model  surface,  the 
projector,  and  the  hydrophone  in  the  tank.  During  any  one  set  of  measurements, 
the  projector  was  held  constant  and  the  hydrophone  moved  in  a circle  in  two 
deg  steps.  Reflection  from  different  areas  of  the  surface  were  obtained  by 
shifting  the  reflecting  surface  horizontally. 

The  projector  was  located  at  a distance  of  1.4  m from  the  surface  measured 
along  the  axis  of  the  projector.  The  hydrophone  was  moved  in  a circle  of  radius 
2.31  m in  such  a way  that  the  hydrophone  and  projector  were  always  pointed  toward 
the  center  of  the  model.  The  mechanical  restraints  in  the  tank  restricted  the 
position  of  the  hydrophone  to  grazing  angles  between  6 deg  and  116  deg. 

*1 

Circular  transducers  whose  active  faces  were  0.125  m and  0.06  m were  used 
to  produce  waves  of  frequency  100  kHz  and  200  kHz,  respectively.  Each  of  these  j 
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transducers  had  a conical  beam  pattern  of  8 deg  at  the  half-power  points.  It 
was  necessary  to  use  a small  hydrophone  to  avoid  spatial  averaging  that  is 
performed  automatically  by  a large  hydrophone.  The  small  omnidirectional 
hydrophone  of  diameter  5.2  mm  used  in  the  experiments  had  a broad  bandwidth 
so  it  could  be  used  at  each  frequency. 


Standard  laboratory  equipment  was  used  to  gate  a continuous  oscillator  for 
signals  of  700  p-sec  duration.  This  produced  an  acoustic  pulse  in  the  water  of 


length  1.05  m which  contained  70  cycles  at  100  kHz.  This  gated  signal  served 
as  the  input  to  a power  amplifier  which  drove  the  projector.  The  received 
signal  was  amplified  by  a tuned  radio-frequency  amplifier  and  displayed  on  a 
cathode  ray  oscilloscope  from  which  the  magnitude  of  the  scattered  pulse  was 
measured.  A representative  photograph  of  the  received  signal  is  shown  on  SLIDE 


' ' 
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The  system  was  calibrated  by  separating  the  projector  and  receiver  a dis- 
tance, rca^,of  3.4  m and  orienting  them  on  axis  toward  each  other.  The  peak-to- 
peak  received  voltage  was  then  V ^ for  this  configuration.  The  transducers 
were  then  arranged  as  shown  in  the  previous  slide.  For  travel  path,  R , 
measured  along  the  acoustic  axes  of  the  transducers  from  the  projector  to  the 
specular  reflection  point  on  the  surface  to  the  receiver,  one  arrives  at  V , the 
peak-to-peak  voltage  of  the  received  signal  in  this  configuration.  The  distances 
and  R are  large  enough  that  one  can  use  the  inverse  square  law  for  pres- 
sure amplitudes  so  long  as  one  is  careful  to  use  the  same  source  signals  for  the 
two  measurements.  One  can  then  calculate  a which  yields  a value  of  unity  when 

fcj 

the  surface  is  a perfect  plane  reflector.  In  the  case  of  a perfect  plane 
reflector,  it  can  be  shown  that  one  would  obtain  a value  of  unity  for  the  scat- 
tering coefficient  if  the  value  J (0,0)  is  normalized  to  a plane  at  right  angles 
to  the  axis  of  the  transducer.  Consequently,  the  quantity  that  should  be  com- 
pared with  the  present  experimental  values  is  a/sintf^. 

Experimental  Results 

Some  of  the  results  of  the  experimental  measurements  will  be  shown  on  the 
following  slides.  On  SLIDES  20  to  25,  the  graphs  show  the  results  of  the  measure- 
ments in  the  plane  normal  to  the  scattering  surface  so  that  0^  = 0 and  b = 0. 
Consequently  Eqs.  1 and  9 can  be  compared  with  this  data.  SLIDE  20  shows  a 
representative  curve  for  ^ = 45  deg  and  a frequency  of  100  kHz.  Measurements 
were  made  and  determined  for  two  deg  increments  of  . Large  excisions  as 
that  occurring  near  55  deg  were  not  unusual.  The  dashed  curve  is  calculated  from 
Eq.  9,  and  it  should  be  emphasized  that  there  are  no  arbitrary  adjustments  of 
parameters  in  this  equation. 
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In  this  slide  and  each  of  the  following  slides,  the  theoretical  quantity 
that  is  plotted  and  labeled  scattering  coefficient  is  not  the  o ^ defined  in 
Eq.  1 nor  defined  in  Eqs.  9 and  10  but  these  quantities  divided  by  (=sin  . 
This  modification  was  justified  in  the  section  on  experimental  techniques.  The 
factor  (1/Cq)  is  a constant  for  each  of  the  graphs  although  it  may  change  from 
graph  to  graph. 

The  position  of  the  bottom  was  moved  relative  to  the  position  of  the  source 


and  receiver, and  seven  more  sets  of  data  similar  to  that  presented  in  this  slide 
were  obtained.  The  average  of  these  eight  curves,  which  gives  a better  statis- 
tical average  for  comparison  with  theory,  is  shown  on  SLIDE  21.  The  curve 
indicated  by  dashes  is  computed  from  Eq.  9,  while  the  curve  plotted  as  circular 
dots  is  computed  from  Eq . 1,  the  formula  derived  from  the  unmodified  theory 
developed  by  Eckart.  This  slide  contains  the  empirical  evidence  mentioned  earlier 
that  justifies  the  boundary  condition  given  in  Eq.  6. 

SLIDE  22  contains  the  experimental  data  and  a theoretical  curve  based  on 
Eq.  6 for  an  incident  grazing  angle  of  90  deg,  i.e.  normal  incidence,  and  a fre- 
quency of  100  kHz.  The  experimental  curve  is  the  average  of  five  measurements  at 
different  positions  on  the  reflecting  surface. 

SLIDES  23  and  24  contain  the  results  of  experiment  and  theory  for  = 40  deg 
and  a frequency  of  100  kHz  and  200  kHz,  respectively.  It  is  pleasing  to  see  that 
Eq.  9 predicts  correctly  the  change  of  cr  with  frequency.  SLIDE  25  contains  the 
results  for  = 30  deg  and  a frequency  of  100  kHz.  The  theoretical  values  based 
Eq.  9 and  plotted  as  a dashed  curve  are  significantly  too  small.  When  it  is 
remembered,  however,  that  the  surface  relief  is  nearly  four  wavelengths,  it  is 
obvious  that  significant  portions  of  the  valleys  are  shadowed  by  the  hills  at 
this  small  grazing  angle.  This  shadowing  produces  a reduction  of  the  topographic 
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relief,  and  one  would  expect  that  the  effective  value  of  h , the  mean  square 

surface  relief,  is  smaller  than  the  measured  value.  This  argument  was  used  to 

2 

justify  recomputing  Eq.  9 with  an  effective  root-mean-square  height  of  0.8(h  ), 

the  value  determined  from  measurements  on  the  model  surface.  This  adjusted 

curve,  which  is  plotted  on  this  slide  with  circular  dots,  agrees  well  with  the 

observations.  This  adjustment  not  only  increases  the  values  of  a , but  it  also 

m 

produces  an  improvement  in  the  shape  of  the  curve. 

In  view  of  the  close  agreement  between  the  theoretical  and  experimental 

curves  seen  in  the  last  several  slides,  one  has  confidence  in  the  predictions 

made  from  Eq.  10  for  other  scattering  angles.  SLIDE  26  contains  a polar  plot 

of  ,9  ) versus  6 for  &.=</>  =40  deg,  f = 100  kHz,  and  the  model  surface 

Mirr  r lr 

relief.  SLIDE  27  is  a similar  plot  of  cl,  versus  6 for  <t>  .=<b  =30  deg.  The 

r M r i r 

2 

effective  value  of  0.8(h  ) for  the  surface  relief  was  used  in  this  calculation 
in  accordance  with  the  agreement  given  above.  These  curves  show  clearly  how 
the  roughness  of  the  surface  broadens  the  scattered  radiation  pattern  in  the 
9^-direction,  and  one  can  see  how,  as  the  grazing  angle  becomes  smaller,  the 
scattered  energy  is  removed  from  the  backward  and  sideward  direction  and  con- 
centrated in  the  forward  direction. 

Experimental  measurements  of  the  scattering  coefficient  as  a function  of 
incident  angle  have  been  made  for  three  other  random  surfaces.  The  experimental 
procedures  used  for  these  surfaces  were  similar  to  those  used  for  the  first 
surface.  The  surface  described  previously  will  be  referred  to  hereafter  as 
SURFACEj.  One  of  the  other  surfaces,  SURFACE  III,  is  similar  to  SURFACE  I 


with  the  exception  that  the  relief  has  been  compressed  by  a factor  of  one-half 


with  respect  to  SURFACE  I.  SURFACE  IV  is  also  similar  to  SURFACE  I except  its 
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relief  has  been  compressed  by  a factor  of  one-fourth  with  respect  to  SURFACE  I.' 
SURFACE  II  has  approximately  the  same  root-mean-square  height  as  SURFACE  I and 
represents  one-fourth  of  SURFACE  I expanded  in  horizontal  dimensions  by  a factor 
of  two.  The  high  frequency  case  as  presented  by  Eckart  and  modified  as  described 
previously  appears  to  be  adequate  for  SURFACE  I when  signals  of  100  kHz  and  200 
kHz  are  used.  For  SURFACE  II,  which  has  the  same  anomaly  height  to  wavelength 
ratio  as  SURFACE  I,  it  is  expected  that  the  high  frequency  limit  will  still  be 
applicable.  For  SURFACE  II  and  especially  SURFACE  III  which  have  root-mean- 
square  heights  of  0.38  cm  and  0.22  cm,  the  applicability  of  the  high  frequency 
limit  is  questionable.  No  comparisons  of  the  high  frequency  limit  or  modifica- 
tions of  this  limiting  case  and  the  experimental  measurements  have  been  made  at 
this  time. 

SLIDE  28  is  a composite  of  the  experimentally  measured  scattering  coefficients 
for  SURFACES  I,  III,  and  IV  for  an  incident  angle  of  40  deg  and  a frequency  of 
100  kHz.  In  the  specular  direction,  it  is  certain  that  the  scattering  coefficient 
for  SURFACE  I,  the  surface  with  the  maximum  relief  of  the  three  surfaces,  is 
below  that  of  the  other  two  surfaces.  However,  from  this  preliminary  data,  no 
definite  statement  can  be  made  on  the  difference  between  the  scattering  coef- 
ficients for  SURFACE  III  and  SURFACE  IV.  SLIDE  29  illustrates  data  similar  to 
that  shown  in  SLIDE  28  with  the  exception  that  the  frequency  is  200  kHz.  Again 
the  scattering  in  the  specular  direction  for  SURFACE  I is  below  that  of  the  other 
two  surfaces,  and  it  appears  that  the  specular  scattering  from  SURFACE  III  is 
slightly  below  that  from  SURFACE  IV.  This  result  might  be  expected,  since  as 
the  surface  becomes  smoother  more  energy  would  be  scattered  in  the  specular 
direction. 

A comparison  of  the  absolute  scattering  coefficients  in  the  specular 
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direction  for  the  three  surfaces  at  100  kHz  and  200  kHz  indicates  little,  if 
any,  variation  with  frequency. 

■ 

SLIDE  30  illustrates  the  measured  scattering  coefficient  for  SURFACE  II, 
the  horizontally  expanded  version  of  SURFACE  I,  for  an  incident  angle  of  30 
deg  and  for  frequencies  of  100  kHz  and  200  kHz.  These  curves  appear  to  indicate 
a definite  frequency  dependence.  SLIDE  31,  however,  does  not  appear  to  indicate 
any  frequency  dependence  when  the  incident  angle  is  40  deg.  Further  analysis 
of  the  data  may  or  may  not  substantiate  this  apparent  frequency  dependence  with 
incident  angle  as  illustrated  by  this  and  the  previous  slide.  An  inspection 
of  the  data  indicates  no  such  frequency  dependence  upon  incident  angle  for  the 
other  surfaces. 

SLIDE  32  illustrates  the  comparison  of  measured  scattering  coefficients  at 
100  kHz  and  for  an  incident  angle  of  40  deg  for  SURFACE  I and  for  SURFACE  II. 

For  receiver  angles  greater  than  30  deg,  these  measured  values  are  approx •mately 
equal.  Theoretically,  in  the  high  frequency  case,  it  would  be  expected  that 
the  scattering  coefficient  for  SURFACE  II  would  be  greater  than  that  for  SURFACE 

[ , • 

I. 

SLIDE  33  illustrates  a comparison  of  the  scattering  coefficient  at  200  kHz 
for  SURFACE  I and  at  100  kHz  for  SURFACE  III.  These  two  situations  involve 
approximately  equal  anomaly  height  to  wavelength  ratios.  For  the  high  frequency 
caae,  ic  would  be  expected  that  the  scattering  coefficient  in  the  specular 
direction  for  SURFACE  III  at  200  kHz  would  be  about  6 dB  greater  than  tn<  sp^c ^lar 
scattering  coefficient  for  SURFACE  I at  100  kHz.  As  is  evident  from  this  slide. 


the  general  trend  for  SURFACE  III  is  somewhat  greater  than  for  SURFACE  I. 
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Conclusions 

Experimentally  measured  scattering  coefficients  as  a function  of  incident 
angle  and  frequency  for  the  rough  surface,  SURFACE  I,  are  in  substantial  agree- 
ment with  predicted  values  using  a modified  version  of  the  short  wavelength 
limit  scattering  formula  developed  by  Eckart.  The  modifications  to  this  theory 
are  related  to  the  specification  of  a boundary  condition  on  the  surface  and  to 
the  technique  of  assuming  root-mean-square  surface  heights  less  than  measured 
values  to  account  for  shadowing  of  portions  of  the  scattering  surface  at  low 
grazing  angles. 

From  the  preliminary  data  available  to  date  for  the  scattering  coefficient 
as  a function  of  incident  angle  and  frequency  for  SURFACES  II,  III,  and  IV,  no 
definite  statements  can  be  made  concerning  comparisons  with  theoretically  pre- 
dicted values.  At  this  time, only  general  trends  and  relative  values  of  the 
scattering  coefficient  in  the  specular  direction  can  be  compared.  Both  the  high 
and  low  frequency  cases  of  Eckart' s theory  will  be  studied  and  possibly  modified 
further  with  subsequent  comparisons  with  measurements. 
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